Electron renormalization of sound interaction with two-level systems in 

superconducting metglasses 
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The crossing of temperature dependencies of sound velocity in the normal and the superconducting 
state of metallic glasses indicates renormalization of the intensity of sound interaction with two-level 
systems (TLS) caused by their coupling with electrons. In this paper we have analyzed different 
approaches to a quantitative description of renormalization using the results of low-temperature 
ultrasonic investigation of Zr4i.2Tii3.sCui2.5NiioBe22.5 amorphous alloy. It is shown that the adia- 
batic renormalization of the coherent tunneling amplitude can explain only part of the whole effect 
observed in the experiment. There exists another mechanism of renormalization affecting only nearly 
symmetric TLS. 
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I. INTRODUCTION 

The well known tunneling model (TM) utilizes only 
two basic parameters for the description of low-tempe- 
rature (T < 1 K) behavior of the velocity v and the 
attenuation T of sound in metallic glasses.EJ A parameter 
Co = pj 2 /pv 2 \p is the density of states of two-level sys- 
tems (TLS), 7 is the deformation potential, p is the mass 
density] defines the scale of variations of v and T in pres- 
ence of TLS. A parameter rj = no^J v 2 k , (rig is the den- 
sity of electron states at the Fermi level, v\ k , is the mean 
square of matrix element of electron-TLS scattering from 
k to k' state) determines the TLS relaxation rate due to 
their interaction with electron environment. According 
to the TM, a TLS contribution to the acoustic character- 
istics is determined by two additive mechanisms - the res- 
onance and the relaxation ones. Under usual experimen- 
tal conditions, uj -C T, where uj is the sound frequency 
(we use the system of units where K = ks = 1) 5 the 
resonance contribution to the variation 5v(T) of sound 
velocity is always negative and represents a straight line 
with the unit slope in coordinates Sv/Cqv versus InT. 
The relaxation contribution is also always negative and 
linear with the slope — 1/2 in the same coordinates. Thus, 
the resulting dependence Sv(\nT) / 'Cqv is expected to be 
a straight line with the slope 1/2. The attenuation of 
sound is determined mainly by the relaxation interaction 
and is virtually independent on T, whereas the resonance 
contribution into Y is small (~ uj/T). 

In superconducting glasses at low enough temperatures 
T <C T c the relaxation interaction is frozen out. This al- 
lows to extract purely resonance contribution and there- 
fore to verify many of TM conclusions. Acoustic m&a-, 
surements, in superconducting metglasses Pd3oZr7 p'El 
Cu30Zr70J.il and (Moi_ x Ru x )o.8Po.2£l carried out more 
than a decade ago, revealed some considerable deviations 
from the predictions of the TM: 

i) the slope of the straight line Sv n (\nT)/CoV in the 



normal (n) phase is about 1/4 whereas the TM canonical 
slope is 1/2; 

ii) at least at high frequencies (HF), the normal state 
line v n (T) crosses the superconducting (s) line v s (T) at 
T cr <C T c . From the TM point of view, this is impossible 
in principle; 

iii) v s (T) is smaller than v n (T) just below T c . This 
effect was observed-in both low frequency (XF) vibrating- 
reed experimentsatm and HF experiments!!! According to 
the original TM, the sound velocity would always increase 
below T c ; 

iv) the sound attenuation reveals an analogous anoma- 
ly: T S (T) exceeds T n (T) within a certain temperature irt. 
terval, which is about T/T c > 0.8 in HF measurements 
and extends down to T/T c > 0.05 in LF experimentsfl 
In contrast, the TM predicts the attenuation to be nearly 
independent on T (with small dT/dT > 0) as long as the 
maximum relaxation rate v exceeds uj. Thus, the atten- 
uation in LF experiments should be insensitive to the 
superconducting transition at all, whereas in HF mea- 
surements r s should either be temperature-independent 
just below T c or decrease in metglasses with low enough 
T c (or at high enough frequencies). 

It was supposed in Refs. |2|,|| that all (or most of all) 
deviations from the TM are related to the electron renor- 
malization of the parameter C with respect to its bare 
value Co- Although possible mechanisms of this renor- 
malization were not discussed in Refs. |]|^, the phe- 
nomenological consideration is rather simple. Indeed, 
assume that C decreases due to the interaction of TLS 
with the electron excitations. As a result, the slope of 
v ra (lnT) decreases also. On the other hand, the bare 
value Co should retrieve far below T c . Therefore, the ra- 
tio of slopes v n (In T) and v s (In T) becomes smaller than 
the canonical TM value 1/2. An additional assumption 
that the parameter C grows more rapidly just below T c 
than the relaxation interaction is frozen out, leads to a 
simple explanation of the items iii) and iv). A connec- 
tion between the item ii) and the renormalization of C 
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is less obvious. Nevertheless, we will demonstrate that 
the crossing between v n {T) and v s (T) at T cr <C T c is the 
most convincing evidence of reduced effective value of C 
in the n-phase in comparison with the s-phase. 

The arguments in favor of the electron renormalization 
hypothesis have been already presented in comparatively 
early theoretical works devoted to both the general prob- 
lem of tunneling with dissipation and a more detailed 
analysis of the TLS interaction with surrounding elec- 
trons (see Refs. HQ and references therein). However, 
any relations allowing to make a comparison with the 
experiment at a quantitative level were not derived in 
these works. 

For the first time, a straightforward theoretical analy- 
sis of the problem of the electron renormalization of the 
sound-TLS interaction in metallic glasses was made in 
Ref. ||. It was argued that one of the reasons for the 
decrease of C in the presence of electrons is an adia- 
batic renormalization of the coherent tunneling ampli- 
tude. Moreover, in order to estimate this effect, it is not 
necessary to introduce any additional parameter since 
the renormalization (Co — C)/Cq — i] 2 /4 is determined 
by the same interaction constant rj (see also Ref. ||). Al- 
though the theorytj gives some opportunity to examine 
its conclusions quantitatively, such procedure was not 
accomplished, probably because of lack of detailed ex- 
perimental data. 

In the present work we test different approaches 
to the quantitative analysis of the sound velocity 
and attenuation in metglasses using experimental re- 
sults obtained on the superconducting amorphous 
Zr4i.2Tii3.8Cui2.5NiioBe22.5 alloy as an example. It is 
shown that the adiabatic renormalization solely does not 
allow to describe all the experimental results, and there 
exists an additional mechanism of renormalization. 



II. EXPERIMENTAL RESULTS 

The alloy under investigation has a high resistance 
with respect to crystallization in the state of overcooled 
melt and remainSp-amorphous at extremely low cooling 
rate (< 10 K s -1 ).E£l This makes it possible to obtain bulk 
homogeneous samples, which suit perfectly the acoustic 
measurements. The .ultrasonic experimental technique is 
described elsewhere .Eil 

Figure 1 shows typical temperature dependence of the 
velocity of the transverse sound wave in Zr4i.2Tii3.s 
Cui2.5NiioBe22.5 in the n- and s-states. The n-state 
measurements were carried out at the magnetic field B = 
1.5-^2.5 T. In accordance with the TM, the curves v(lnT) 
represent almost straight lines in both the n- and the 
deep s-state. The growth of v s below T c reflects freezing 
out of the relaxation component and agrees with the TM 
conception, with the constant Co = (2.85 ± 0.05) • 10~ D 
determined from the slope of v s (lnT) at T < 0.3 K. 
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FIG. 1. Temperature variations of the velocity of trans- 
verse sound in Zr4i.2Tii3.8Cui2.5NiioBe22.5 alloy in the su- 
perconducting and normal phases. Inset: normalized velocity 
of transverse (upper trace, C„ = 6.94- 10~ 6 ) and longitudinal 
mode (lower trace, C„ = 2.75 • 10 -6 ) near T c . The curves 
were aligned in the normal phase. 
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FIG. 2. Temperature variations of magnetic susceptibil- 
ity and sound velocity in Zr4i.2Tii3.sCui2.5NiioBe22.5 in the 
vicinity of T c . Open and solid circles: B = and B — 1.5 
T. Solid lines: calculations for rj = 0.65, £& = 1.2, Ub = 0.5, 
R s = 0.14, T c = 0.83 K. Thick vertical mark shows noise level. 
Experimental data was smoothed by adjacent averaging. 

There are also obvious deviations from the TM: the 
ratio of slopes in the n- and the s-phases differs from its 
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canonical value 1/2 and is close to 1/4, and the curves 
v(T) for both phases intersect at some temperature T cr . 
Such effects have been observed before in Pd3oZr7 .B The 
value of T cr is frequency-dependent; particularly, T cr (62 
MHz) fa 0.055 K and T cr (186 MHz) ps 0.11 K were found. 

The inset to Fig. 1 shows the variations v(T) for trans- 
verse (t) and longitudinal (1) sound normalized over cor- 
respondent slopes in the n-phase. Obviously, these two 
dependencies virtually coincide, whereas absolute veloc- 
ity variations are sufficiently different. 

Figure 2 shows temperature dependencies of v s and v n 
in the vicinity of T c and the diamagnetic response of the 
sample on the ac magnetic field H = 10~ 6 T at the fre- 
quency of 22 Hz. The magnetic susceptibility x was mea- 
sured simultaneously with v s that provides coincidence of 
temperature scales for both measurements. The presence 
of two steps in x(T) indicates that the sample contains 
at least two phases with different temperatures of the su- 
perconducting transition: T cml pa 0.9 K and T cm2 pa 1,0 
K. An increase of the ac field H up to 10~ 5 T leads to 
complete suppression of the anomaly in the diamagnetic 
response at T cm2 , although the jump at T cml survives up 
to H = 10~ 4 T. 
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FIG. 3. Normalized attenuation of the transverse sound 
wave versus temperature in Zr41.2Ti13.sCu12.5Ni10Be22.5- 
Open and solid circles: B — and B — 2.5 T. Solid lines: 
calculated dependencies for r\ = 0.65, T c = 0.83 K. Inset: to 
determination of the superconducting energy gap and the pa- 
rameter r). The results are presented for three values of the 
background level of attenuation: —0.3%, 0% and +0.3% of 
the total change of attenuation between the n- and s-phase, 
from bottom to top curve. 

It is of interest to note that the temperature T c pa 0.83 
K, at which one can first register a nonzero difference 
v s — v n , coincides neither with T cm i nor with T C m2. This 



was interpreted in Ref. O, as a fingerprint of possible 
gapless superconductivityH within the temperature re- 
gion between T c and T cm i . It is known that the magnetic 
scattering, which is the most prevalent reason of the gap- 
less regime, also reduces the energy gap when the latter 
opens. However, we found the energy gap in our alloy to 
be close to the BCS value (see below) that allows us to 
reject this interpretation. Apparently, the diamagnetic 
anomalies are related to some surface phases with higher 
transition temperatures. 

The temperature dependence of the sound attenuation 
is shown in Fig. 3. The experimental data are normalized 
over T n (T c ) which can be determined from the variations 
of sound amplitude between T c and deep s-state. The be- 
havior of r(T) /r(T c ) does not show any noticeable differ- 
ence with similar dependencies in other superconducting 
amorphous alloys and reflects evolution of the relaxation 
contribution to the attenuation, in accordance with the 
TM conception. 



III. QUALITATIVE CONSIDERATION 

Before making quantitative estimations, we shall dis- 
cuss qualitatively a possible origin of the peculiarities of 
the sound velocity in superconducting metglasses. 

It is naturally to associate the crossing of v n (T) and 
v s (T) with a growth of C in the s-phase as a result of 
suppression of the electron renormalization at energies 
smaller than the superconducting gap. To validate this 
assumption, we address the expression for-the resonant 
contribution of TLS to the sound velocityO 



P 



CEtanh(E/2T) 



E 2 



dE. 



(1) 



In the simple case of energy- independent C — Co, in 
order to avoid a formal logarithmic divergence of Eq. (1) 
at the upper limit, the velocity variations are usually 
considered with respect to some arbitrary reference tem- 
perature Tq\ 



v(T) - v(T ) 
v(T ) 



-Coin—, T>lo. 
Jo 



(2) 



However, if the value of C varies with energy and/or 
temperature, Eq. (2) is inapplicable even for qualitative 
estimates since the reference value v(Tq) may also change 
with C. To account correctly for the changes of C, it is 
necessary to analyze the complete integral of Eq. (1) by 
introducing a cutoff energy E m (say, of the order of melt- 
ing temperature or glass transition temperature). In the 
case of C = const, Eq. (1) can be approximated within 
the logarithmic accuracy by the following piecewise- linear 
dependence (line 1 in Fig. 4): 



5v(T) 
Cv 



ln(u/E m ), 
HT/E m ), 



T < uj 

T> LJ 



(3) 
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Here we neglect insignificant small variations in v at T < 
lo: a quadratic fall near T — and a shallow minimum 
at w_= 2.2T arising from the analytical solution of Eq. 
(1)0 
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FIG. 4. Schematic representation of temperature depen- 
dence of sound velocity in superconducting glass. Lines 1, la 
and lb correspond to (Sv(T) / Cov) Yes in the superconducting 
phase. Line 2 represents (6v(T)/Cov) lea in the normal phase. 
Line 3 is the total (8v(T) /Cqv) in the normal phase. See text 
for details. 

Since the parameter Co is determined for TLS in the 
absence of electron excitations, the low-temperature part 
of line 1 in Fig. 4 depicts the resonance contribution to 
the sound velocity at temperatures much lower than T c . 
Let us assume that the resonance interaction in the n- 
phase is determined by a renormalized constant C < Co 
for all energies below a certain threshold E^. By making 
use of a simple step-like approximation (C = Co(l — R), 
< R = const < 1 for E < Ek and R = in the oppo- 
site case) , one can plot the resonance contribution in the 
n-phase by line 2 which is located over the correspond- 
ing line for the s-phase. Thus, the resonance interaction 
increases the sound velocity in the n-phase at T < Ek if 
C decreases and vice versa. This basic conclusion cannot 
be derived from Eq. (2), where "the dc component" is 
lost. 

The relaxation contribution to the sound velocity is 
much less sensitive to the renormalization effects and be- 
comes negligibly small at T <j w/r] 2 . This allows us to 
write down the following approximate expression: 



(HT)\ = 1 
V C v Ad 2 



0, 

\n{u}/Trf 



T < oj/rf 1 
T > uj/t] 2 



(4) 



The total TLS contribution to the sound velocity in 
the n-phase is schematically shown in Fig. 4 with a line 
3 which was obtained in the following way. At T < Ek 
the line 3 contains a piece of line 3n drawn with a slope 
1/2 - R through the point T = ui/rj 2 at the line 2. At 
T > Ek, the line 3 is a piece of line 3s drawn with a slope 
1/2 from the point T = oj/tj 2 at the line 1. A transition 
range arises at temperatures T sa Ek- 

Let us now discuss the evolution of v{T) at the super- 
conducting transition. First, we consider a low-tempera- 
ture range T < 0.3 K. If E k < 2A S (0) (here and below 
A S (T) is the superconducting energy gap), the renormal- 
ization of C should be frozen out. Thus, the resonance 
contribution can be depicted by a low-temperature part 
of the straight line la. It is clear that the intersection 
of lines 3 and la is possible only under the condition 
C < Co {R 7^ 0)- From the geometry of Fig. 4 it is easy 
to estimate T cr : 

( 1- S) In T cr = R\nE k + (l-R-S) ln(c^/n 2 )—A(T CT ), (5) 

where S is the resulting slope of v n (In T) . The parameter 
A is introduced to account for possible shift of a back- 
ground level of the sound velocity in the s-phase with 
respect to the n-phase normalized over Co and will be 
discussed below. It can be seen from Eq. (5) that T cr 
grows along with Ek, to, and with the decrease in 77. 

At Ek > 2A S (0) the renormalization of C is only par- 
tially frozen out for E < 2A S (0). Therefore, the sound 
velocity in the deep s-state can be depicted as a part 
of the straight line lb (see Fig. 4) . The estimate of T cr 
by Eq. (5) is also valid in this case, only Ek should be 
replaced by 2A S (0). 

Along this line of reasoning, we can also qualitatively 
explain the behavior of v and r at the superconduct- 
ing transition. Below T c , the electron renormalization 
rapidly reduces, and the effective C grows providing the 
decrease in v and the increase in T. However, a com- 
petitive effect arises simultaneously: the rate v of the 
TLS relaxation on electrons falls and therefore changes 
v and r in the opposite direction. Thus, if the phonon 
relaxation predominates near T c , the effective v changes 
weakly, and the sound velocity will decrease (correspond- 
ingly, the attenuation will increase) below T c , as it was 
observed before-CTB If the electron relaxation prevails (for 
materials with lower T c like our system), the changes of 
v and r near T c may have any sign, depending upon the 
relations between T c , Ek and u. 

In principle, one can propose an alternative explana- 
tion of the crossing. In our previous consideration, wc 
silently assumed the count level of the TLS contribution 
into the sound velocity to be the same for both the n- and 
s-states. Generally, this is not the case, and the sound 
velocity changes at the superconducting transition with 
no account of the TLS-related mechanisms. For exam- 
ple, in pure metals a decrease in electron viscosity below 
T c leads to the change of dislocation contribution to the 
sound velocity of the order of 10 -5 (Ref. [l6]) which is com- 
parable with the TLS contribution but is undoubtedly 
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absent in an amorphous metal. A more general mech- 
anism is the change of electron contribution into elastic 
moduli of a metal in the s-phase. In disordered metals 
with a short electron mean free path, this change is usu- 
ally small (~ 10~ 6 ) but in certain cases, for instance, in 
A-15 compounds close to structural instability, it may 
achieve much larger values ~ 10~ 4 (Ref. ^) of arbitrary 
sign. If we accept such a scale for the decrease in the 
electron contribution in the s-state of our sample, the 
crossing will arise without any renormalization effects. 
In its turn, the anomalous slope ratio may be attributed 
to enhanced density of states of asymmetric TLS play- 
ing the principal role in the sound attenuation. Although 
the latter assumption contradicts the basic TM postulate 
about constancy of p within a wide range of tunnel pa- 
rameters, we can not reject straight away the discussed 
alternative without additional argumentation presented 
below. _ 

Thermodynamic treatmentLj shows that the electron 
contribution variations in the s-state are independent on 
the sound frequency and lead to a jump in derivatives 
dvi/dT at T = T c for both longitudinal and transverse 
modes proportional to d 2 T c /de 2 where are corre- 
sponding deformations. A small jump of v\ itself can 
be also expected at T = T c . As the temperature de- 
creases, the electron contribution changes as the density 
of a superfiuid condensate, so that its variations become 
negligibly small at T <C T c . By making use of Eq. (5) 
at R = and the measured values of S ~ 0.28, Co, 
r\ f» 0.65 and T cr (62 MHz), the corresponding shift of the 
sound velocity between the n- and the deep s-state can 
be estimated as Sv/v ~ 5 • 10~ 5 . This value is compa- 
rable with the resulting velocity change 3 • 10~ 5 for t- 
mode (see Fig. 1) between T c and the maximum in v(T), 
i.e., the electron and the TLS contributions appear to 
be of the same order but have opposite signs. However, 
since the normalized TLS contribution is independent 
on the polarization, the data presented in the inset in 
Fig. 1 demonstrate that the same independence should 
take place for the electron contribution, i.e., the condi- 
tion {l/j 2 )d 2 T c /de 2 « (l/-f 2 )d 2 T c /de 2 must be satis- 
fied. The latter does not follow from theory and can be 
only a result of random coincidence that is hardly pos- 
sible. Thus, we conclude that the scale of temperature 
variations of the TLS contribution much exceeds that 
connected with the electron mechanisms. Furthermore, 
since the velocity shift A is frequency-independent, the 
following expression derived from Eq. (5) 

(l-S)ln^p\ = (l-R-S)ln^ (6) 

shows that in the absence of renormalization (R = 0) T cr 
must be proportional to ui that contradicts our experi- 
mental data. In that way, the absence of such propor- 
tionality is the most clear evidence of renormalization of 
the parameter C irrespectively of whether a certain addi- 
tional sound velocity shift between n- and s-phases exists 
or not. 



IV. SOME RESULTS OF THE TUNNELING 
MODEL 

In this Section we present a brief overview of basic 
results of the TM, which describe the behavior of v(T) 
and r(T) in glasses with account of the dependence of 
C on tunnel parameters and were used in our numerical 
calculations. We also discuss modifications introduced 
into given relationaibr a more exact account for the TLS- 
electron coupling.m] 

The main postulate of the TM is a statement of the ex- 
istence of double- well potentials in glasses with the tunnel 
coupling between wells. The density of states of TLS is 
constant in the space of parameters £, In Ao where £ is 
the asymmetry of the double- well potential and Ao is an 
amplitude of the coherent tunneling. In order to deter- 
mine the response of the TLS ensemble on an external 
field, it is necessary to perform an averaging over £ and 
InAo, which appears to be more convenient in variables 
E = y/g 2 + Aq and u = Aq/E. In this representation, 
the TLS density of states is independent on E: 



g(E,u) 



P 



= g(u). 



(7) 



The relationships which determine the TLS contribu- 
tion to the sound velocity and attenuation reaca: 



^v(T) 



tanh 



(1) 



C(e, u)g(u)u 2 du 



(8) 



6v(T) 



i r Em/T ds r 1 



rcl 



2 7o cosh'(e/2)7 



C(e,u)g(u) 



x (1 - u 2 )- 



-,du, 



(9) 



Tv 



rcl 



E m /T 



de 



coslr (e/2) 



C(e,u)g(u) 



x (1 -u 2 ) 



LU 2 + U J 



;du. 



(10) 



In Eqs. (8)-(10) we introduced e = E/T. 

In general case, the relaxation rate v is determined by 
both the electrons and phonons, but for T < 1 K the 
phonon contribution can be neglected. In the original 
TM the TLS-electran interaction is considered within a 
perturbation theoryEJ over the parameter ?/ 2 , which does 
not affect the splitting of energy levels, and all specific 
features of the metglass, in comparison with amorphous 
dielectrics, reduce only to the appearance of a new relax- 
ation channel having the rate 
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v = ^f-u 2 TJ(s). 



(11) 



In the n-phase J(e) = J n {e) = (e/2) coth(e/2), v 
r] 2 Tu 2 , and the relaxation interaction is essential for ,&L1 
T > to. In the s-state it is necessary to use a functionta 



e'(e' - e) - A 2 
y/{e' - e) 2 - A 2 



6(( e '- e) 2 - A 2 )sign(e'- e) + (e — e) , (12) 



J{e'~e) 
f{-e) 



where f(x) is the Fermi function, Q(x) is the step Heavi- 
side function, and A = A S (T)/T. This integral coincides 
with J n (e) for e 3> 2 A, it has a discontinuity at e = 2 A, 
and J s (e,A) -> 2/(A) for e < 2A. A rapid fall in J s 
below T c leads to freezing out of the relaxation interac- 
tion when the maximal relaxation rate (u = 1) becomes 
smaller than u). 

A more complicated picture was revealed beyond the 
perturbation theory.Hu Just at T = in the n-phase, the 
bare amplitude of the coherent tunneling Ao is renormal- 
ized due to an adiabatic part of the interaction between 
the TLS and electrons: 



A* 



UJQ 



(13) 



where loq is of the order of the Debye energy. 

In the n-phase for T ^ an ensemble of TLS can be 
divided into three energy intervalsH 

1. T <E* = v/£ 2 + A* 2 - coherently tunneling TLS. 

2. E* < T < AE/ttt] 2 - incoherently tunneling TLS 
having the tunneling amplitude 



Acx A (2^r/wo) ?72/4 



(14) 



and the energy splitting of E = y £ 2 + A 2 . In new vari- 
ables E* , E and u*, u within the intervals 1, 2, respec- 
tively, the relations Eqs. (7)-(10) hold. 

3. T > 4E/irr) 2 - low -energy TLS are relevant. In this 
region, the amplitude of incoherent tunneling is also A 



of Eq. (14) but the factor (1 



in Eqs. (9) and (10) 



is absent because the incoherent transitions between the 
broadened levels take place with energy variations even 
in the symmetric case. The corresponding relaxation fre- 
quency varies as 



2 

7T?7 2 



V% OC ttTu — — 

- 2 J(e) 



(15) 



One can think that a part of TLS with E < \JojT 
should decrease its contribution to T and (Sv/v) le i due 
to the fall of v% for small e. However, a numerical anal- 
ysis shows that this fall is compensated by the growth 



of influence of the symmetric TLS. As a result, partial 
contribution to T and (8v/v) Ie \ from the interval 3 vir- 
tually does not change in comparison with the original 
TM. The contribution of TLS from the second interval re- 
mains the same also. Only the contribution from the co- 
herently tunneling TLS, undergoing the adiabatic renor- 
malization, experiences an essential change. The main 
postulate of the TM concerning the constancy of p in the 
space of variables £, In Ao remains valid. However, g(u*), 
along with the parameter C, acquires an additional fac- 
tor (1 — 77 2 /4) under transformation to the variables E* , 
u* because of a nonlinear dependence between Aq and 
Ao, Eq. (13). The latter changes rapidly at the super- 
conducting transition to the linear one of the type of Eq. 
(14), where As(T) substitutes for T, and all TM relations 
are restored. u 



V. ANALYSIS OF THE EXPERIMENTAL DATA 

A. Determination of A s (0) and r\ 

For the frequencies used in our experiments, a rapid 
freezing of relaxation interaction begins at the tempera- 
ture well below T c (see Fig. 3). In this case the renormal- 
ization of C is also frozen out and the sound attenuation 
T s is described by Eq. (10) with C = Co. According to 
Eqs. (10) and (11), the low-temperature part of T S (T) 
should be a straight line in coordinates \n(T s (T)/T) ver- 
sus T" 1 : 



T fl (T) _ 2irrf 

r„(T c ) 



3w 



Te 



-A.(0)/T 



(16) 



This allows us to use the sound attenuation for a sim- 
ple evaluation of A s (0) and rj from its low-temperature 
dependence plotted in the inset to Fig. 3. Since this con- 
struction is very sensitive to the reference level of the at- 
tenuation, we also present two additional curves for level 
variations of ±0.3% of the whole signal change between 
the n- and s-states. Within this range, a large enough 
temperature interval exists where each curve can be well 
approximated by a straight line with the slope determin- 
ing A s (0) = 1.45 ±0.05 K. If we accept T c = 0.83 K, this 
value agrees well with the BCS relation A s (0)/T c = 1.76 
which was used in all further calculations. 

The value of r\ = 0.55 ± 0.15, determined by cross- 
ing of the approximating straight lines with the ordinate 
axis for different curves in the inset to Fig. 3, reveals a 
large spread due to its exponential dependence on the 
position of the crossing point. A more accurate estimate 
of rj can be obtained from the numerical analysis of the 
attenuation within the whole temperature region of Fig. 
3. By matching the most sharp part of T S (T), calculated 
from Eq. (10) at T c — 0.83 K, with the experimental de- 
pendence, we found ?y = 0.65 ± 0.05, in agreement with 
previous rough estimate. 
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B. Sound attenuation near T c 

The analysis performed above shows that the tem- 
perature dependence of the sound attenuation in 
Zr4i.2Ti 13 8 Cui2.5Ni 10 Be 2 2.5 can be rather well described 
by the original TM. However, within some temperature 
range just below T c , the behavior of T s reveals anomalies 
which find no explanation in the original TM. 
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FIG. 5. Comparison of sound attenuation near T c at 54 
MHz with model calculations. Open and solid circles: 5 = 
and B — 1.5 T. Lines: related calculation for the original 
TM (set 1 of lines) and the TM with account for the adi- 
abatic renormalization of the coherent tunneling amplitude, 
rj — 0.65, Eb = 1.2 (set 2). Thick vertical mark shows noise 
level, data was smoothed by adjacent averaging. 

Figure 5 shows variation of T(T) in the vicinity of T c 
obtained with higher resolution than in Fig. 3. The main 
peculiarities of r s (T) can be discovered by means of a 
comparison of the experimental curves with those calcu- 
lated in the frame of the original TM (line set 1 in Fig. 
5). According to the calculation, a fall in the attenua- 
tion begins just at T c with growing slope at low temper- 
atures. The experimental dependence has quite different 
behavior: r s (T) does not vary at T c within the exper- 
imental resolution and tends to exceed r„(T) at lower 
temperatures. A more prominent excess of T S (T) over 
T n (T) starting just at T c has been previously observed 
in PdsoZryo alloy.fj 

The proximity of A s (0)/T c to the BCS value, as well as 
the lowered temperature of superconducting transition of 
the main volume of the sample in comparison with T cm \ 
(see Fig. 2) shows that the explanation of T S (T) features 
by the magnetic depairinglia is irrelevant. Let us now 
discuss the applicability of the concept of the electron 
renormalization of C to the description of the behavior 



ofr s (T). 

Indeed, Fig. 1 shows that the renormalization of C re- 
ally takes place that follows from the crossing of v s (T) 
and v n (T) at T cr <C T c . However, the ratio of slopes of 
w(lnT) in the n- and s-states points out that the scale of 
renormalization is rather large: 5C/Co ~ 0.22, i.e., more 
than twice possible maximal contribution rj 2 /A ~ 0.1 of 
the adiabatic renormalization of Ao . Hence there should 
be an additional mechanism of renormalization. More- 
over, an incomparability between the scale of SC/Cq and 
the magnitude of the anomalies of T s (T) indicates that 
this mechanism involves TLS which do not contribute to 
the relaxation attenuation at T ~ T c . Recall that the 
main contribution to the attenuation comes from TLS 
with u opt « lu or u opt ~ \f^h 2 T < 1. 

One of possible additional mechanisms of the renormal- 
ization is associated with the fluctuational modulation of 
a barrier in the double-well potential. This mechanism 
affects only almost symmetric TLS (u ~ 1)0 and does 
not give a contribution into T S (T). 

According to Ref. ^, the adiabatic renormalization of 
A should involve all TLS with E >T. This condition 
places coherently tunneling TLS into the range where a 
cutting factor in the denominator of Eq. (10) is relevant. 
In spite of that, their partial contribution to T(T) can be 
essential in the scale of Fig. 5. 

In the original TM, the value of T(T C ) for lu < T c is 
close to 1/2 for the scale used in Fig. 5. The decrease in C 
shifts r n (T) towards smaller values. In order to analyze 
the shift of the experimental dependence with respect 
to the calculated one, we need the accuracy better than 
1% for the absolute value of the attenuation, which is 
beyond the accuracy of our experiments. Therefore, in 
Fig. 5 we discuss only the relative position of T n (T) and 
T S (T) curves, which was measured much more precisely. 

For numerical calculations we used the following en- 
ergy dependence of the parameter C: 

S- =l-^Q(es b ) [l + (2/(A)-l)e(2A- e )]-B. (17) 

This approximation is quite reasonable, since for e* fi. 
1 the coherent amplitude AJ decreases exponentiallya. 
Here e b ~ 1 is a fitting parameter which confines the 
range of coherently tunneling TLS. The last factor in the 
second term takes into account that for e < 2 A the con- 
tribution comes only from the normal excitations. The 
last term in Eq. (17) takes into consideration an addi- 
tional renormalization due to the symmetric TLS; its ori- 
gin will be discussed below. A contribution of B to the 
sound attenuation is negligibly small and at this stage we 
assume B = for the sake of simplicity. 

The results of simulation are plotted in Fig. 5 (lines 
2). The difference between lines In and 2n reflects a 
contribution from the adiabatic renormalization of Ao to 
the sound attenuation in the n-state. We matched al- 
most temperature- independent part of T S (T) at rj = 0.65 
with measured curves and obtain quite reasonable value 
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of Eb = 1.2 ± 0.1. The calculated dependence of T S (T) 
varies similarly to the predictions of the original TM just 
below T c . Then T S (T) undergoes a break, changing the 
sign of dT/dT at T = 2A s (T)/e b . These features arise 
due to exploiting a step approximation in Eq. (17). The 
superconductivity does not affect the renormalization of 
C if 2A S (T) is smaller than the value of E = Te b . Obvi- 
ously, a smoothed energy dependence of the cutoff factor 
in Eq. (17) would decrease the variation of T S (T) at T c 
and eliminate the break. The same result would also oc- 
cur due to a possible broadening of the superconducting 
transition in an amorphous sample. 

In this way, the evolution of T S (T) in the vicinity of 
T c is determined by two factors: a fall due to the de- 
crease in the relaxation rate v and a growth because of 
freezing out of the C renormalization. The first factor 
is frequency-dependent in contrast to the second one. 
Therefore, the resulting variations of T s should also de- 
pend on frequency. When u) decreases, the temperature 
range where T s (T) > F„ (T) should be extended and vice 
versa. In particular, if ry, T c , and e b are fixed, the in- 
crease in frequency by an order of magnitude (see, for 
example, Ref. || where the measurement frequency was 
about of 600 MHz) has to mask the action of the second 
factor utterly. This frequency increase leads to T s (T) al- 
ways lower than T n (T). Besides, dT S (T) / dT grows as the 
temperature decreases. However, the experiments in Ref. 
^ were carried out in metglass with T c — 2.6 K where v is 
essentially determined by phonons and depends weakly 
on the state of electron subsystem. Under these circum- 
stances, the freezing out of the renormalization should 
give even stronger effect than observed in our case. 

C. Sound velocity 

The resonance contribution to v(T) is completely de- 
termined by coherently tunneling TLS with s > 1. 
Therefore, the account of adiabatic renormalization of 
Ao in the n-state will lead to (Sv/CQv) Tes = (1 — 
T] /4)lnT independently on the magnitude of e b in Eq. 
(17). From Eqs. (9) and (17), we estimate a magni- 
tude of the relaxation contribution as (5v/Cqv) tc \ — 
-1/2 [l - ?7 2 /4(l - tanh(e 6 /2))] InT. Using the values 
of i] and Sb obtained before, we get (Sv/Cov) n — 0.42 InT 
for the total change in the n-state, whereas the slope 
of the experimental dependence (Sv/Cov) n — 0.28 InT 
differs from the original TM coefficient 0.5 much more. 
Thus, the mechanism of the adiabatic renormalization of 
Ao can solely explain less than a half of the whole effect, 
and, as it was mentioned above, an additional origin of 
the electron renormalization of C has to exist. It must 
affect mainly the symmetric TLS which do not partici- 
pate in the relaxation attenuation. Besides, both mecha- 
nisms can be considered as additive because the scale of 
C renormalization is small. 

The authors of Ref. ^ studied the effect of electron 



density fluctuations on the barrier height of the inter- 
well potential. It was argued that below some critical 
temperature almost symmetric TLS and a surround- 
ing electron cloud can form a strongly-correlated (bound) 
state similar to the Kondo state. This effect has an en- 
ergy threshold, E < Ek(Tk). Fluctuations also lead to 
the renormalization of the tunnel amplitude similar to 
Eq. (13): 

A^ = A (T/£>f\ (18) 

where D is of the order of the Fermi energy. The expo- 
nent m depends on rj and is close to 0.1-^0.2 for r\ ~ 0.65. 

The spectrum transformation of Eq. (18) does not re- 
sult in the renormalization of the parameter C since it 
does not change the density of the TLS states within the 
space of new variables u, E. These questions were not 
considered in Ref. |7[ Nevertheless, one can conclude that 
the renormalizing factor would depend onu= Aq/.E in 
Eq. (18) since, according to Ref. |7j, A — > A if u — > 0. 
This nonlinear relation between Ao and Ao means the 
effective renormalization of TLS density of states similar 
to the adiabatic renormalization. One can also expect a 
reduced value of the deformation potential in the bound 
state. 

Thus, the slope (Sv/Cov) ICS is determined by the re- 
lation C/Co = 1 — rj 2 /A — R s where R s describes the 
contribution from symmetric TLS in n-phase. The latter 
does not change the slope of (Sv/Cov) re i, therefore R s is 
not a fitting parameter. Its value for given rj is unam- 
biguously determined by the resulting slope S (with the 
account for a small correction related to the contribution 
of adiabatic renormalization of Ao in (Sv/Cav) rc i). Par- 
ticularly, the estimates presented above give R s ~ 0.14. 

For computing, we model B in Eq. (17) by simplest 
step function, introducing a conventional boundary u b of 
the "symmetric" TLS: 

B = _^_e( u - tt5 )[i + (2/(A)-l)e(2A-e)]. (19) 

The meaning of two last factors in Eq. (19) is clear 
from preceding discussion. The first factor represents the 
"real" renormalization by the symmetric TLS because 
the renormalizing correction appears with the weight 
a^/I — under the integration over u in Eq. (8). 

A comparison between calculated and experimental de- 
pendencies for Ub = 0.5 is presented in Fig. 6. Here we 
use the following procedure. The experimental points 
for the frequency of 62 MHz were taken from Fig. 1 and 
normalized by Co . It is impossible to measure mutual po- 
sition of v(T) for different frequencies with necessary ac- 
curacy of 10~ 7 . However, according to Eq. (3), the value 
(Sv(T)/v) Tes does not depend on frequency for uj <C T. 
That is why the experimental points in the s-state for 
both frequencies 186 MHz and 62 MHz are placed to- 
gether at temperatures far below T c . The positions of 
the calculated dependencies are determined by the upper 
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limit of integration in Eq. (7). Finally, one can also com- 
pare them with experimental dependencies after match- 
ing of (Sv(T)/Cv) ICS in the s-state at T < T c . 




0,1 1 
T(K) 

FIG. 6. Comparison of the experimental change of sound 
velocity for different frequencies and magnetic fields (circles: 
62 MHz, triangles: 186 MHz; open symbols: B = 0, solid 
symbols: B — 2.5 T) with the calculations. The following set 
of parameters is used: rj = 0.65, Et = 1.2, = 0.5, R 3 = 0.14, 
T c = 0.83 K. 

Figure 6 shows satisfactory agreement between the cal- 
culated and the experimental dependencies. Some differ- 
ence arises only in the temperature range from T c to 0.4 
K (see also Fig. 2). The calculated curve in the s-phase 
is noticeably steeper at T ~ T c than the experimental 
one. An estimate shows that a more smooth energy de- 
pendence of the adiabatic renormalization of Ao leads 
only to insignificant decrease in the slope of (5v/Cqv) s 
at T ~ T c . The most probable reason for deviations of 
the computed dependencies from the experimental ones 
near T c is the smearing of superconducting transition and 
a contribution from_small thermodynamic corrections to 
v(T) in the s-phase.EZI A small contribution can also come 
from a residual influence of phonon relaxation which may 
cause the excess of the calculated slope of (5v/Cov) n over 
the experimental one near T c (Fig. 2). 

The account for the symmetric TLS shifts up the re- 
sulting dependence (Sv/Cov) n at R s In 1.36u( without 
changing the slope of (5v/Cqv) ic i, as follows from com- 
putation of Eq. (9) with renormalization of Eq. (19) for 
oj/rfT <C u\ <C 1. This shift plays the same role as 
the parameter A in Eq. (5). Therefore, according to Eq. 
(5), the fitting of Ub is reduced in fact to matching of 
the computed T cr with the experimental value at given 
rj, R S) and S. The fitted value of Ub — 0.5 is smaller 
than Ub ~ 0.7 accepted in Ref. R. However, the latter 



is determined by the bare magnitude of Ao , whereas the 
integration in Eq. (8) is performed over the renormalized 
variable u. Therefore, according to Eq. (18), Ub will be 
smaller than the bare value. 

The account for the thermodynamic correction also 
will lead to an increase in Ub\ formally, Ub can be made 
as much close to 1 as one likes. Unfortunately, there is 
no clear way to separate the electron contribution on the 
background of the TLS effects. Note that a weak temper- 
ature dependence in Eq. (18) will introduce a correction 
into the slope (8v(lnT)/Cov) rc \ too. Then the value R s 
in Eq. (19) should be decreased by a factor of (1 + to). 

Thus, our assumption about an additional mechanism 
renormalizing the contribution of almost symmetric TLS 
to the sound velocity gives a satisfactory description of 
the sound velocity behavior in both the n- and s-states 
with a help of only one extra parameter Ub- Note that 
a thorough calculation of T S (T) for Ub = 0.5 with the 
additional parameter B in Eq. (16) leads to practically 
the same dependencies as shown in Fig. 5 for B = 0. 

VI. CONCLUDING REMARKS 

Using the results of the acoustic measurements ob- 
tained on the superconducting glass Zr4i.2Tii3.gCui2.s 
NiioBe22.5 we have carried out a quantitative analysis of 
different theoretical approaches to the electron renormal- 
ization of sound interaction with TLS. A convincing evi- 
dence for the renormalization is the crossing of the lines 
v n (T) and v s (T) at some temperature T cr <C T c , in com- 
bination with the absence of proportionality between T cr 
and uj. Assuming the simplest model renormalization of 
the interaction parameter C in the space of tunnel vari- 
ables, it is possible to describe quantitatively the behav- 
ior of the sound velocity and attenuation exploiting the 
original tunnel model. It is sufficient to use the adiabatic 
renormalization of the coherent tunneling amplitudetffl to 
fit the dependence Y S {T) with the experiment. However, 
the behavior of the sound velocity can be described only 
with the help of an additional mechanism of the renor- 
malization affecting only almost symmetric TLS. This 
is the main result of our consideration. The additional 
mechanism can be presumably related to-irebuilding of 
the interwell potential due to fluctuations^ but this ap- 
proach has not been developed enough to consider it to 
be incontrovertible. 

The analysis carried out in this work allows us to eval- 
uate several parameters using the experimental depen- 
dencies of T(T) and v(T). We would like to emphasize 
that most of them are not fitting parameters in a com- 
mon sense because it is not necessary to vary the whole 
set of them simultaneously to determine each of them. 
We utilize the following sequence of parameter evalua- 
tion. First of all, from the temperature dependence of 
the sound attenuation in the s-state we determine the 
energy gap A s (0) and the parameter rj. Using the lat- 
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ter together with the experimental T S (T), v n (\nT), and 
T cr , we evaluate the parameters Sb, R s , and conse- 
quently. In order to evaluate any parameter, we utilize 
only the values which have been already found during 
previous steps. In this way, rj is the central parameter 
which determines all others. Therefore a question may 
rise: could an error in the determination of rj lead to es- 
sential redistribution of the contributions from different 
mechanisms? Under assumption that the electron contri- 
bution into v s (T) does not change below T c , the needed 
value of Ub ~ 0.3 is nonrealistically small already for 
r\ = 0.67 -j- 0.68. On the other hand, for r\ < 0.6 the best 
fit yields lower T c . In this way, the studied alloy does not 
give a large choice for the variation of the parameters. 

In principle, it would be possible to manage without 
the concept of additional influence of electrons on the 
symmetric TLS in the analysis of elastic properties of our 
metglass and reduce all effects to the adiabatic renormal- 
ization of the tunnel amplitude, if we accept X] w 1 assum- 
ing its preliminary estimate from low-temperature "tail" 
of the attenuation (Sec. V, A) to be unreliable. How- 
ever, in this case the overall shape of calculated T S (T) 
can be adjusted to the experimental dependence only at 
T c = 0.9 K, revealing noticeable deviations from the ex- 
perimental data at T <C T c . Besides, the experimental 
value of T cr could be obtained only under assumption 
that in s-state at T — T cr the electron contribution is 
lowered by Sv/v ss 1.5 • 10 -5 (this estimate follows from 
Eq. (5)) that is comparable with the TLS contribution. 
Within this approach, the resulting variations of v(T) 
near T c would be determined by three equipollent mech- 
anisms: the change of the TLS relaxation rate, the freez- 
ing of the adiabatic renormalization, and the evolution 
of electron contribution. In our opinion, it is impossible 
to expect proper mutual compensation of their partial 
contributions, for both the longitudinal and transverse 
modes, providing regular variation of v(T) observed at 
the superconducting transition. 

It is interesting to notice that the parameter Ek which 
appears in the fluctuation modefl is not so important in 
our evaluation process. Our analysis does not demand 
also the absence of this parameter. Most probably, the 
condition Ek > 2A S (0) is satisfied in the alloy used in our 
investigation because, as follows from Eq. (5), the intro- 
duction of Ek < 2 A s (0) at constant T cr should be accom- 
panied by a decrease in Ub- So we believe that there is no 
need to use Ek for the description of low-temperature fea- 
tures of the sound velocity and attenuation in our alloy. 
The introduction of Ek can be probably more fruitful for 
the analysis of the metglass elastic properties at higher 
temperatures. 
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